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Abstract. Two different deterministic microscopic traffic flow models, which 
are in the context of the Kerner's there-phase traffic theory, are introduced. In 
an acceleration time delay model (ATD-model), different time delays in driver 
acceleration associated with driver behaviour in various local driving situations 
are explicitly incorporated into the model. Vehicle acceleration depends on local 
traffic situation, i.e., whether a driver is within the free flow, or synchronized 
flow, or else wide moving jam traffic phase. In a speed adaptation model (SA- 
model), vehicle speed adaptation occurs in synchronized flow depending on driving 
conditions. It is found that the ATD- and SA-models show spatiotemporal 
congested traffic patterns that are adequate with empirical results. In the ATD- 
and SA-models, the onset of congestion in free flow at a freeway bottleneck 
is associated with a first-order phase transition from free flow to synchronized 
flow; moving jams emerge spontaneously in synchronized flow only. Differences 
between the ATD- and SA-models are studied. A comparison of the ATD- and 
SA-models with stochastic models in the context of three phase traffic theory is 
made. A critical discussion of earlier traffic flow theories and models based on the 
fundamental diagram approach is presented. 
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1. Introduction 

Theoretical studies of freeway traffic flow dynamics is one of the rapid developing fields 
of statistical and nonlinear physics (see the reviews P El S El Ej j the book j2| , and 
the conference proceedings |51l^ llU[[mil2lll3[ll4lll5j '). For a mathematical description 
of freeway traffic flow, a huge number of different microscopic and macroscopic traffic 
flow models have been introduced. In macroscopic models, individual dynamic vehicle 
behaviour is averaged, i.e., these models describe dynamics of average traffic flow 
characteristics like average vehicle speed and density (see e.g., [TCI HTl IT51 [T^ I^ U. 
Microscopic traffic flow models describe individual dynamic vehicle behaviour, which 
should simulate empirical spatiotemporal features of phase transitions and congested 
patterns in freeway traffic. In this article, we restrict a consideration of microscopic 
traffic flow models only. 

There are two types of microscopic traffic flow models: Deterministic models and 
stochastic models PElElSElElIZi- In deterministic models, some dynamic rules 
of vehicle motion in traffic flow are responsible for spatiotemporal features of traffic 
patterns that the models exhibit. Contrastingly, stochastic models, in addition to 
dynamic rules of vehicle motion, exhibit model fluctuations, which play a fundamental 
role for traffic pattern features. 

There are at least two classes of deterministic traffic fiow models El El El El- In 
the first class, driver time delays in vehicle acceleration (deceleration) a are explicitly 
taken into account. An example is the classic model of Herman, MontroU, Potts, 
and Rothery l23|: If the vehicle speed v, or the speed difference between the vehicle 
speed and the speed of the preceding vehicle ve, or else the net distance g (space gap) 
between vehicles changes, then the driver accelerates (decelerates) with a time delay 

r m 

a{t + T)^f{v{t),ve{t),g{t)). (1) 

Based on Gazis, Herman, and Rothery (24j have developed a microscopic traffic 
flow model, which is capable of describing traffic beyond of instabilities; steady state 
solutions of this model lie on a one-dimensional curve in the flow-density plane (the 
fundamental diagram) (see the review by Nagel et al. [HI for more detail). Recall 
that steady state solutions are hypothetical model solutions in which all vehicles move 
at the same time-independent speed and the same space gap between vehicles. One 
of the mathematical descriptions of this model class first proposed by Nagatani and 
Nakanishi |2SI and further developed by Lubashevky et al. [2E1 reads as follows 

da f{v{t),ve{t),g{t)) - a{t) 
dt T 

In both models |24l 1251 ESj , steady state model solutions in the flow-density plane lie 
on the fundamental diagram. 

There is also another class of deterministic microscopic models in which the vehicle 
speed satisfies the equation P El El El El ^ 

^ = ^iv{t),v,it),g{t)). (3) 

Examples are optimal velocity (OV) models of Newell 1271, Whitham |2HI, Bando, 
Sugiyama et al. jJHl, and the intelligent driver model (IDM) of Treiber and 

X It should be noted that transferring the information delivered from one vehicle interacting with 
the neighbour ones requires to deal carefully with a complex averaging process by derivation of a 
macroscopic traffic flow model. The related mathematical theory is developed in Ref. I21II22I . 
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Helbing PU] E]- Steady state solutions of this model class that obviously satisfy 
the conditions (f)(v,vi,g) = and v = V£ lie on the fundamental diagram in the 
flow-density plane. 

If functions and model parameters in the models Q and ||3J) are chosen in an 
appropriated way, then there is a range of vehicle density in which steady state 
model solutions for free flow are unstable. This instability, which should explain 
the onset of congestion, leads to wide moving jam emergence in free flow (F^J 
transition) 01311 0^1. 

However, as explained in the book [7j, the above models that are in the context 
of the fundamental diagram approach, as well as all other traffic flow models reviewed 
in m 121 El El m cannot explain the fundamental empirical feature of traffic 
breakdown, i.e., that the onset of congestion in free flow at a bottleneck is associated 
with a local first-order phase transition from free flow to synchronized flow (F^S 
transition) |321 133L rather than with an F-^J transition. For this reason, Kerner 
introduced a three-phase traffic theory. In this theory, there are three traffic phases: 
(i) free flow, (ii) synchronized flow, and (iii) wide moving jam. 

The first microscopic models in the context of three-phase traffic theory 
introduced in 2002 are stochastic models jSSlESI- As in empirical observations |33l I34j . 
in these models wide moving jams emerge spontaneously only in synchronized flow 
(S— *J transition), i.e., the models exhibit the sequence of F^S^J transitions leading 
to wide moving jam emergence in free flow; in addition, the models show all 
types of congested patterns found in empirical observations |2H1 IMl I2Z1 OHl Q 12^1 ■ 
Recently, some new microscopic models based on three-phase traffic theory have been 
developed j4()l 1411 1^ . However, there are no deterministic models in the context of 
three-phase traffic theory, which can exhibit the F^S— >J transitions found in empirical 
observations and the diagram of congested patterns of three-phase traffic theory |Mir7|. 
In stochastic models 1231211 (23 EH Hi driver time delays in acceleration (deceleration) 
are simulated mainly through the use of model fluctuations. Therefore, a development 
of deterministic models based on three-phase traffic theory is important for a more 
realistic theory of car following behaviour. 

In this paper, two deterministic microscopic three-phase traffic models are 
presented. In an acceleration time delay model (ATD-model for short; Sect.jJll, an 
explicit description of driver time delays in vehicle acceleration (deceleration) is used. 
In a speed adaptation model (SA-model for short; Sect.|21l, vehicle speed adaptation 
occurs in synchronized flow depending on driving conditions. In Sects. 01 and El 
we show that these models exhibit the F— >S^J transitions and congested patterns 
associated with results of empirical observations. In addition, a stochastic SA-model 
is introduced and compared with the deterministic SA-model of Sect. El In Sect. El the 
deterministic microscopic three-phase traffic models of Sects. (21 and El are compared 
with earlier deterministic models and a critical discussion of models in the context of 
the fundamental diagram approach is performed. 

2. Acceleration Time Delay Model 

2.1. Driver Behavioural Assumptions and Empirical Basis of ATD- Model 

A deterministic three-phase traffic flow model with driver time delays (ATD-model) is 
based on the following empirical features of phase transitions and congested patterns 
as well as driver behavioural assumptions of three-phase traffic theory (Sects. 2.3, 2.4, 
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and 8.6 of the book |7|): 

(i) In synchronized flow, a driver accepts a range of different hypothetical steady 
states with various space gaps g at the same vehicle speed v, i.e., steady states of 
synchronized flow cover a two-dimensional region in the flow-density plane. 

(ii) To avoid collisions, in the steady states a driver does not accept the vehicle 
speed that is higher than some safe speed (denoted by Vs{g, vi)) that depends on the 
speed of the preceding vehicle ve. In contrast with earlier models in which a safe speed 
determines a multitude of steady states on the fundamental diagram PI31EIH1E1E3|, in 
the ATD-model the safe speed determines the upper boundary of the two-dimensional 
region for the steady states in the flow-density plane Jj . 

(iii) If a driver cannot pass the preceding vehicle, then the driver tends to adjust 
the speed to the preceding vehicle within a synchronization gap G{v, vi), i.e., at 

g < G{v, V,) (4) 

a speed adaptation effect occurs. The synchronization gap determines the lower 
boundary of the two-dimensional region for the steady states in the flow-density 
plane. In the ATD-model, the speed adaptation effect is modelled through a driver 
acceleration K{v, ve){v — ve) adjusting the speed to the preceding vehicle under the 
conditions K(v, ve) is a sensitivity. 

(iv) In traffic flow with greater space gaps, a driver searches for the opportunity 
to accelerate and to pass. This leads to driver over-acceleration, which is modelled 
through a driver acceleration A{V^^'^'"^'> (g) ~ v) adjusting the vehicle speed at 

g > G{v, ve) (5) 

to a gap-dependent optimal speed in free flow V^^'^°°'>{g), where A is a sensitivity of 
this effect. 

A competition between the speed adaptation effect and driver over-acceleration 
simulates a first-order F^S transition leading to the onset of congestion in real traffic 
flow (see explanations in Sect. 2.4 in (Tj). 

(v) In empirical observations, due to an F— >S transition there is a maximum 
point of free flow associated with the maximum density pmax"*, maximum flow rate 
Qmax°\ and maximum speed w^j'n"^ given by the formula v^^^^"^^ = qma-^^ / Pma^'' (Sect. 
2.3 in [7j). This maximum point is modelled through F— i-S transition, which occurs 
already due to infinitesimal local perturbations in steady states of free flow associated 
with the optimal speed in free flow V^^^'^''\g) at the density pmax''. 

(vi) In high density flow, a driver decelerates stronger than it is required to 
avoid collisions if the preceding vehicle begins to decelerate unexpectedly (driver 
over-deceleration). In the ATD-model, the over-deceleration effect, which explains 
and simulates moving jam emergence in synchronized flow, is modelled by a driver 
time delay r in reduction of a current driver deceleration (denoted by t = t['^°'^\v)). 
The longer t['^'^'^\ the stronger the over-deceleration effect. In empirical observations, 
the lower the synchronized flow speed, the greater the probability for moving jams 
emergence (Sect. 2.4 in [7]). For this reason, t['^'^'^\v) is chosen to be longer at lower 
speeds than at higher ones. 

(vii) At the downstream front of a wide moving jam or a synchronized flow region, 
a driver within the jam or the synchronized flow region does not accelerate before the 
preceding vehicle has begun to accelerate. In the ATD-model, this effect is modelled 
through the use of a mean driver time delay in acceleration at the downstream front 
of the synchronized flow region, which depends on a time delay in driver acceleration 
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(denoted by r = t^'^ ) and on the sensitivity K{v, vg) at w < vg. At the downstream 
front of a wide moving jam, a mean time delay in acceleration from a standstill v — 
within the jam should be longer than the mean driver time delay in synchronized 
flow [?]• To simulate this longer mean time delay in vehicle acceleration, in addition 
with two mentioned above model effects, a vehicle within the jam does not accelerate 
before the condition 

9 > .g^r^ ^ (6) 

is satisfied, in which gma™'' is the maximum space gap within the wide moving jam 
phase. 

(viii) Moving in synchronized flow of lower speeds, a driver comes closer to the 
preceding vehicle than the synchronization gap G. In empirical observations, this 
self-compression of synchronized flow is called the pinch effect (Sect. 12.2 in |7|). 
In the ATD-model, the pinch effect is simulated through the use of two model 
assumptions. Firstly, a time delay in reduction of a current driver acceleration 
(denoted by t = q-^'^'^'^') increases if the speed decreases. Secondly, the sensitivity 
K{v, Vg), which describes the speed adaptation effect (item (iii)), is chosen a.t v > vg 
different from K{v, vg) at v < vg (item (vii)). Specifically, K{v, vg) at v > vg is 
chosen to be smaller at low speeds than at higher ones. As a result, at lower speeds 
vehicles choose smaller space gaps than the synchronization gap G. 

(ix) At the upstream front of a wide moving jam or a synchronized flow region, a 
driver begins to decelerate after a time delay denoted by r = rjf'^^''. This delay time 
should describe realistic velocities of deceleration fronts in congested traffic patterns. 



2.2. Main Equations 



An ATD-model reads as follows: 
dx 

Hi 

dv 

'dt 

da 
'dt 



a, 




at g > G and g > .g^l" ■* , 
at g < G and g > 5^ax\ 
at < g < gidax ^ , 



(7) 
(8) 

(9) 



where x is the vehicle space co-ordinate; g = xg — x — d; the lower index £ marks 
variables related to the preceding vehicle; all vehicles have the same length d, which 
includes the minimum space gap between vehicles within a wide moving jam; a^^'"^''\ 
a^'^^"-' , and a(j™) are vehicle accelerations (deceleration) in the free flow, synchronized 
flow, and wide moving jam phases, respectively. If the condition |(SJ| is satisfied, then 
a vehicle moves in accordance with the rules for free flow. Within synchronized flow 
associated with the condition gniax' < <7 < G{v, ve), the vehicle tends to adapt the 
speed to the preceding vehicle. Within a wide moving jam, the space gap is small, 
specifically g < (?ma'x\ and the vehicle decelerates. § 

§ Since the vehicle speed v cannot be negative, the following condition is also used for Eqs. JSJ, 1^ : 

a{t) > at v{t) = 0. (10) 

To satisfy this condition in numerical simulation, the acceleration a{t) is replaced by the value 
max(a(t), 0) if v{t) = at time t. 
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2.3. Driver Acceleration 

The accelerations (decelerations) a(^'^°°\ a*^''^'^-' , and a*^-''''"^ are found from the condition 

a(phasc) ^^j^^(^^^(g(phasc)^ ^^^^^ 

the superscript "phase" in means either "free" , or "syn" , or else "jam" for the 
related traffic phase; flg is a deceleration related to safety requirements; amin and flmax 
(omin < 0, flmax > 0) are respectively the minimum and maximum accelerations for 
cases in which there are no safety restrictions. In Hll(l . functions and 
^(jam) associated with driver acceleration within the related traffic phase ~ free flow, 
or synchronized flow, or else wide moving jam - are determined as follows: || 

V, ^;,)^A(F(f-<=)(5)-z;) + 

K{v, ve){ve-v), (12) 

V, V,) - Amin(Vl7,")(g) -v,0) + 
K(v, vi){vi - v), (13) 

a(jam)(^) ^ _^(jam)^^ (14) 

Here Knax"' (5) is a gap-dependent maximum vehicle speed in synchronized flow; i^O''™) 
is a sensitivity. 

2.4. Safety Conditions 

Safety deceleration with a deceleration can be applied, if the vehicle speed becomes 
higher than the safe speed Vs{g, ve). We use safety deceleration found from the 
condition: 

as{g, V, vg) = As{vs{g, ve) - v), (15) 

where Ag is the sensitivity related to safety requirements. 

The speed Vs{g, vg) in (|15|l is found based on the safety condition of Gipps |43j : 

v,T, + vl/{2b,) <g + vj/{2bs), (16) 

where Tg is a safety time gap, 6s is a constant deceleration.^ We use an approximated 
formula for Vs{g, vi) derived from (|16|l in [Appendix A| which enables us to write 
as (5, V, ve) fTB|l as follows 

as(g, V, vi ) = A^^\vi){glTg - !>) + Kg(vt)(vi - v), (18) 

where 

Ai^Hve) = AgTg{Tg + ve/{2h)y\ (19) 
Kg{vi) = Asin + ve/{2h)KTs + v,/{2hg))-\ (20) 

To is a constant. 

II In the article, large enough flow rates on the main road are considered at which congested patterns 
can occur at a bottleneck. For this reason, in 1121 K is chosen to be independent on g in the free flow 
phase. At considerably smaller flow rates in free flow, specifically, if g increases, K in 03 should 
tend towards zero when g ^ G. 

f Note that Eqs. of the ATD-model can also be written without the term a^^™' as follows 

da^j (a(f-°) -a)/r at g > G, , 
dt \ (a(=y") -a)/T at g < G. ^ ' 

In 1171 . the speed Vs{g, vi) in 1151 is found based on the Gipps-condition III6I when g > (/majT' ^^nd 
the speed Vs{g, v{) = when g < gmasT'- In the latter case, the formula 1151 with Vs(g, vi) = plays 
the role of vehicle deceleration within the wide moving jam phase. 
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2. 5. Physics of Driver Time Delays 

In Eqs. the time delay r is chosen as 

Ts atas<min(0, niax(a(P'^'''''=\ flmm), a), 
f otherwise. 

Here, Tg is a short driver time delay associated with a finite driver reaction time that 
must be taken into account in the cases when the driver should decelerate unexpectedly 
to avoid collisions; f is a time delay in other traffic situations, which is chosen different 
depending on whether the vehicle accelerates or decelerates: 

r^'^'^'^) at a > 0, 



r(d<=c) at a < 0. 



In turn, t^^^^^i and r^'^'^'^' in (|22|l depend on the acceleration a 

'i'^'^''^ at a < a(P'^'*'"= 
otherwise. 



-^^^^^ = <i V) "^r "":" (23) 



{^('^°'^') of ^ ^ ^(phaso) 
^|dcc) otherwise. 

The driver time delays Tq^'^'^^ , Tq'^'^'^^ , t['^°'^^ ^ and r}'*'^'^-' in (|23() , H24() are associated 
with human expectation of local driving conditions, in particular, with spatial and 
temporal anticipation of a driver in accordance with local adaptation to those traffic 
situations in which the driver takes into account both the current and expected future 
behaviour of many vehicles ahead (see also Sect. 12. 

^^acc) mean time delay when a driver starts to accelerate or wants to increase 

the acceleration. This can often occur at the downstream front of a wide moving jam 
or a synchronized flow region, i.e., when the speed in traffic flow downstream of the 
vehicle is higher than the current vehicle speed. In these cases, after the preceding 
vehicle has begun to accelerate, the driver also begins to accelerate, however, after a 
time delay to have a desired time gap to the preceding vehicle. 

Tg'^'^'^^ is the mean time delay when the driver starts to decelerate or wants to 
decelerate harder in cases in which the driver approaches a region of a lower speed 
downstream. 

^(acc) (corresponds to situations in which the driver accelerates currently but wants 
either to stop the acceleration or to reduce it. Thus, t^^'^'^'^ is the mean driver time 
delay in interruption or reduction of driver acceleration in cases in which the driver 
recognizes that current acceleration is greater than a desired acceleration in the current 
driving situation. 

^(dcc) (corresponds to situations in which the driver decelerates currently but wants 

(dec) 

either to stop the deceleration or to reduce it. Thus, Ti is the mean time delay in 
interruption or reduction of driver deceleration in cases in which the driver recognizes 
that current deceleration is more negative than a desired deceleration in the current 
driving situation. 

2. 6. Model of Road with On-Ramp Bottleneck 

"Open" boundary conditions are applied on the main road of the length Lq. At the 
beginning of the road free flow conditions are generated for each vehicle one after 
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another at equal time intervals Tin = 1 /im where q-m is the flow rate in the incoming 
boundary flow. To satisfy safety conditions, a new vehicle appears only if the distance 
from the beginning of the road {x = xy,) to the position x = xi of the farthest upstream 
vehicle in the lane exceeds the distance veTg + d. The speed v and coordinate a; of a 
new vehicle are v = ve and x = xi^, respectively. After a vehicle has reached the end 
of the road, it is removed; before this, the farthest downstream vehicle maintains its 
speed. In the initial state {t = 0), all vehicles have the same initial speed v = l/^^''"') (g) 
and space gap g, and qin = v/{g + d). 

An on-ramp bottleneck on the main road is considered. The on-ramp consists of 
two parts: (i) The merging region of the length Lm that begins at a: = Xon- Within 
this region, vehicles can merge onto the main road from the on-ramp. (ii) The part 
of the on-ramp lane of length L,- upstream of the merging region at which vehicles 
move according to the model equations for a homogeneous road with the maximum 
speed Wfioc, on = 90 km/h. At the beginning of the on-ramp lane the flow rate to the 
on-ramp gon is given as the flow rate on the main road qi^. 

The following rules are applied for vehicle merging within the merging region. A 
speed V is calculated corresponding to formula 



Here superscripts g'^ and g are space gaps to the preceding vehicle and the trailing 
vehicle on the main road, respectively; v'^ and v~ are speeds of the preceding vehicle 
and the trailing vehicle, respectively; 7, g^™'' and Aw^^' are constants, where g'^^^ is 
the minimum gap at which vehicle merging is possible, Av^^^ describes the maximum 
possible increase in speed after vehicle merging. Note that the finite increase Awr^^ in 
the vehicle speed (|25|) is used to simulate a complex driver behaviour during merging 
onto the main road, especially in synchronized flow: In some cases, before merging 
the driver has to accelerate abruptly, to adjust the speed to the speed of the preceding 
vehicle. 

If the conditions 12611 are satisfied, then the vehicle merges onto the main road. 
After merging the vehicle speed v is set to ii (|25|1 and the vehicle coordinate does not 
change. If the conditions H26|) are not satisfied, the vehicle does not merge onto the 
main road. In this case, the vehicle moves in the on-ramp lane until it comes to a stop 
at the end of the merging region. 

2.7. Model Functions and Parameters 

Model functions and parameters are shown in Tables ^ and respectively. As 
explained in Sect. 12.11 driver time delays r}'^°'^' and t['^'^'^^ are chosen to be functions 
of the vehicle speed; additionally, the synchronization gap G{v, vi) and sensitivity 
K{v, Ve) are chosen to be asymmetric speed functions depending on whether the 
vehicle speed v is higher or lower than the speed ve. Explanations of the function 
K{v, Ve) have been made in item (vii) and (viii) of Sect. [TTI 

Speed dependence and an asymmetric function for the synchronization gap 
G{v, Ve) are explained by driver behaviour as follows. The synchronization gap is 
the space gap at which a driver adapts its speed to the speed of the preceding vehicle. 
Firstly, the synchronization gap is an increasing function of speed: The lower the 




(25) 



(26) 
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Table 1. ATD-model functions 



Synchronization gap 



G{v, ve)^vmax{0, r(">"^)(z;) + Vi){v - ve)) 
r ^(acc) nt V <ve 

r(syn)(^) ^ r^('^y")(l _ 0,S5{v/Vo)^) 

Sensitivities 

, f K'-'""'^ a.tv <vi 
^ 1 if (d-) at « > 

A(z;) = (l + exp((i;/t-c-l)/e))-i 

Characteristic speed functions 
F(f>--)(g) = T/(.9) 

= Votanh((5 + 2)/(For)) 

Speed dependensies of time delays 

(dec)/ ^ _ J 0.5 S at W > Wc 



(acc)/ N f 0.57 s at w > 



0.7 s otherwise 
0.57 s at V > Vc 
0.87 s otherwise 



Table 2. ATD-model parameters 



Vq = 33.3 m/s (120 km/h), T = 0.9 s, 

A = 0.5 s-\ if(^'='=) = 0.8 s-i, = 1 s~\ 

y^(syn) ^ s, = 0.5 sVm, Kf'i^") = 0.55 sVm, 

gt"^^ = 0.95 m, r'^^^' = 1 s, r^'^ = 0.75 s, n = 0.4 s, 
Omax = 1 ni/s^, a,„i„ = -1 m/s^, 
As = 1.25 s"\ 6s = 2 m/s2, = 1 s, To = 0.42 s, 
Lo = 25 km, Xb = — 5 km, Xon = 16 km, 

= 300 m, L,. = 500 m, 7 ^ 0.22, Av^^^ = 8 m/s, ff^;;,^ ^ 



speed, the smaller the maximum gap at which the driver can comfortably move in 
synchronized flow. Secondly, if w < ve, the driver accelerates and he/she can start 
speed adaptation at a smaller space gap than in the opposite case v > vi. The 
function T^^'y'^^{v) is used to have a difference in vehicle space gap in steady states of 
free flow and synchronized flow at a given flow rate. This space gap difference, which 
is used for simulation of a first-order F— ^S transition, tends towards zero when the 
density in free flow approaches the maximum point for free flow pmax ^ (figures ^ (a) 
and (c)); see also item (v) of Sect. l2.1|l . 
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Figure 1. Model steady states for the ATD-models: (a) - In the space-gap— 
speed plane, (b) - In the flow— density plane, (c) - Steady states and the line J 
(explanation of the line J see in the book 7 ). (d) - ATD-model with separated 
steady states in free flow and synchronized flow. 



2.8. Steady States 

In steady states, all vehicles have the same speed v = vg and the same space gap g, 
and all accelerations and their time derivatives are zero, and the density p and the 
flow rate q are related to the space gap g and to the speed v by the obvious conditions 

p^l/{xe-x) ^1/ig + d), pv ^ v/{g + d). (27) 

According to and formulae for (g), V^^Tx"' {g) (Table[TJ) for steady states, 
we get 

v = V{g) at5>GH and5>g(j:-), (28) 

v<V{g) at.g<G(t;) and5>g(f-), (29) 

^; = at g<gt^\ (30) 

V < vs{g, v). (31) 

According to (|28|I - H31|I . the model steady states consist of the curve v = V{g) 
H28|l a,t g > 5^;°°^ (curve F in figure ^ (a)) for free flow, a two-dimensional region in 
the space-gap-speed plane for synchronized flow determined by inequalities in H29() . 
((33, and the line v = at g < gmax^ (PH]) for wide moving jams (figure H (a)). 5^-°°^ 
is the minimum space gap in free flow found as a solution of the set of the equations 
V = V{g) and g = G{v, ve) at v = ve. 

The two-dimensional region for steady states of synchronized flow is limited by the 
following boundaries: the boundary U, the curve L, the curve v — V{g) at g < gmin°^ 
and the horizontal line g = g^^\ The boundary U is associated with the safe speed, 
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i.e., this is determined by the condition H31|l when it is an equahty. This leads to the 
condition for the boundary U 

g = vT,. (32) 

The boundary L is found from the condition that the vehicle space gap is equal to the 
synchronization gap 

g = Giv). (33) 

In the flow-density plane, free flow (curve F in figure ^(b)) is found from 

q^pVFip) (34) 

at p < pma°x^ wheTeVpip) = V{g) \g=p-i-d, Pma°x^ = (5'mir''+'^) A wide moving jam 
is associated with the horizontal line q — at p[^i^^^ ^ P ^ Pmax (figure ^(b)), where 

Pmt™^ — (ffmax^ + c?)~^, Pmax = . Thc boundaries of a two-dimensional region for 
steady states of synchronized flow are: the upper line U determined by the condition 
g = (1 — pd)/Ts, the lower curve L determined by the condition pG{q/ p) = 1 — pd, the 
curve lIH^ at p> pmax\ and the vertical line p — p'"^^'^^ 



3. Speed Adaptation Model 

3.1. Empirical F—^S^J Transitions as Physical Basis of Speed Adaptation Model 

The fundamental hypothesis of three-phase traffic theory, which postulates that 
hypothetical steady states of synchronized flow cover a two-dimensional region in the 
flow-density plane, is also one of the basic hypotheses of the ATD-model presented 
above (figuresn(a) and (b)). In contrast with the ATD-model, in a speed adaptation 
model (speed adaptation model, SA-model for short) hypothetical steady states of 
synchronized flow are associated with a curve (curve S in figures [21 (a) and (b)), 
i.e., they cover a one-dimensional region in the flow-density plane. The curve S is 
associated with an averaging of an infinite number of steady states of synchronized fiow 
to one synchronized flow speed for each vehicle space gap. A gap dependence of the 
average speed in synchronized flow steady states on the curve 5* is denoted by Vav^"^ (g) . 
The basis hypothesis of the SA-model is associated with the sequence of F^S— >J 
transitions, which determine moving jam emergence in empirical observations [331 17] . 

Note that as in the models and theories in the context of the fundamental diagram 
approach in the SA-model steady state model solutions cover a one- 

dimension region(s) in the flow-density plane. However, in the models and theories 
reviewed which claim to show spontaneous moving jam emergence, the 

F— >J transition governs the onset of congestion. This is inconsequent with empirical 
results 13 • In contrast, in the SA-model the onset of congestion is associated with an 
F— >S transition, whereas moving jams occur spontaneously only in synchronized flow, 
in accordance with empirical results. 

The SA-model is simpler than the ATD-model. However, due to this simplification 
the SA-model cannot show some features of congested patterns of the ATD-model 

We have also studied another version of the ATD-model in which there is a separation of steady 
states in free flow and synchronized flow in the flow-density plane, i.e., the maximum point for 
free flow Pm^x' is related to the intersection point of the line U and the curve F (figure (d)). 
Simulations of this version of the ATD-model show qualitatively the same features of phase transitions 
and congested patterns as those discussed in Sect. HI 
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speed [km/h] density [vehicles/km] 




Figure 2. Model steady states for SA-models: (a, b) - In the space-gap— speed (a) 
and the flow-density plane (b) for the SA-model 11351 - 11411 . (c, d) - The SA-model 
l35i - l4UI . I42i . (e, f) - Other possible SA-models (see |Appendix B) . Dashed parts 
of curves for model steady states in (a-f) are associated with unstable model 
steady states, (g, h) - Steady states and the line J for the SA-model (a, b) and 
the SA-model (c, d). 
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(Sect. EH), which are observed in empirical observations. The purpose of the SA-model 
is to simulate an F^S transition and features of the sequence of F^S^J transitions, 
as observed in empirical observations [331 [7], in a simple way. This confirms an 
assumption of three-phase traffic theory that if rather than the fundamental hypothesis 
the hypothesis about the F— >S^J transitions is the basis of a mathematical model, 
then the model can show and predict some important empirical features of the phase 
transitions (see footnote 4 of Sect. 4.3.4 in [J]). 

In the SA-models, an F— >S transition is modelled through two effects: (i) 
Discontinuouty of steady speed solutions (figures|21(a), (c), and (e)) or their instability 
(curve FS in figure |21 (f)) in the vicinity of the maximum point of free flow I'^j™''^ 

Pmax''. (ii) The speed adaptation effect is modelled through the term K{v, ve){v — Vi) 
that adjusts the speed to the preceding vehicle in synchronized flow. 

Moving jam emergence is simulated through an instability of some of the 
synchronized flow model steady states associated with the curve 14v'*'"^(g). This 
instability occurs in synchronized flow at lower speeds and greater densities (i.e., 
smaller space gaps). The associated critical density and speed of the synchronized 
flow steady states are denoted by pcr''^ and Wcr''', respectively (figurc[2Jl. To simulate 
this instability, as in the ATD-modcl (item (viii) of Sect. I2.l|l . in the SA-models the 
sensitivity K(v, vg) at v > vi is a, decreasing speed function. Similarly with the ATD- 
model, to simulate the mean time delay in acceleration at the downstream jam front 
in the SA-model, a vehicle within the jam does not accelerate before © is satisfied 
(item (vii) of Sect. ED. 



3.2. Basic Equations 

There can be different possibilities for a separation of steady states of free flow and 
synchronized flow in SA-models, which all exhibit qualitatively the same features of the 
F^S^J transitions. To illustrate this, here we consider two variants of SA-models; 
in |Appcndix B| other possible variants of SA-models are discussed. All these variants of 
the SA-models exhibit very similar features of phase transitions and spatiotemporal 
congested traffic patterns that are associated with the same physics of these SA- 
models. 

A formulation for the SA-model reads as follows 

'^'^ V, (35) 



dt 



dv 
'dt 



^(trcc) at w > Wjjji^ and g > g, 



max 



a(«y") at < z;^r^ and 5 > si^^x \ (36) 
aO™) atO<g<g^Al^\ 
The vehicle acceleration a — dv / dt in (|36|) is supposed to be limited by the maximum 
acceleration Omax, i-e., in H3()|) 

Here and below the associated designations of functions and parameters have the same 
meaning as those in the ATD-model (Sect.^)). 



3.3. Vehicle Acceleration 



Functions a!'^'^°°\g, v, vg), a^'^^'^'>{g, v, vg), and a^^^™\v) in 137|l are chosen as follows 
a(f'='=)(ff, V, Vi) = A(f"^)(T/(*'-°°)(5) -v) + 
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K{v, ve){ve-v), (38) 
a(^>'")(g, V, V,) = A^^y'^^V^^y'^^g) - v) + 

K{v, ve){ve-v), (39) 

a(jam)(^,) ^ _^(jam)^^ (40) 

Two versions of functions Vav'^"'' (g) in (|39|1 that lead to two different versions of 
the SA-models are considered: 

V^r\g) = ~g{9)/Tt''\ (41) 

and 



(42) 



where g(g) — g — (?ma'x Tav''"', Vi and c are constants. 
3.4- Steady States and Model Parameters 

In the SA-models, in accordance with H3()|l there are three isolated curves for steady 
states of the SA-models associated with the three traffic phases: free flow, synchronized 
flow, and wide moving jam (figures [3 (a) and (b)). 

Steady states of free flow are related to a curve v — Vf(p) and formula 134|l (the 
curve F in figures |21 (a)-(d)) associated with the condition 

z; = y(f™)(5) at«>z;^™). (43) 

Steady states of synchronized flow are related to a curve S in the space-gap-speed 
plane (figures |21 (a) and (c)) given by the condition 

v = V^y^\g) at«<^;^:r^ and5>g(C^'. (44) 
In terms of the flow rate q and density p, the formula H44|l reads 

q=pvsip) at pir ^ < p < (45) 

where Vs{p) = Vav'^"^(5) |g=p-i_d, Pmfn'' = (ffmax'' + rf)"\ Smax' IS found from the 
equation Kil^'^"-' (g^^S^ ) = v'^^n^ ■ 

In the case of the function T4v'^"''(5) given by (|41() . the formula 145() yields the 
equation for a curve S with a negative slope in the flow-density plane (figure [5] (b)) 

9 - (1 - P/p!i!n^')/ri7"' at p(:f.r' < p < fJ^. (46) 

When the function Vav^"'' (g) is given by formula (|42|l , the curve S has a maximum in 
the flow-density plane (flgure|21 (d)). 

Steady states for a wide moving jam are the same as those in the ATD-model, 
i.e., they are given by a horizontal line 

9 = piZ^ < p < p^,^ (47) 

in the flow-density plane (figures |2(b) and (d)). 

Parameters of the SA-models are shown in Table El 
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Table 3. SA-model parameters 

V{g) = Vot&nh{g/{VoT)), Vq = 33.3 m/s (120 km/h), T = 0.85 s, 
a„,ax = 2 ms-2; = 0.4 

^(syn) ^ 0.1 S-\ = 2.2 S-\ 

K{v, Vi) is given in Tabled 

/i:(dcc),^^<) is given in Table [Hat Vc = 10 m/s, e = 0.07, 

5min ~ ymax , 7 — U.ZO. 



,-1 



In the SA model with function V^7''\g) EU, 

we use 

= 22.22 m/s (80 km/h), ifj^^^' = 0.95 s-^K^^'^ = 0.64 s" 
^(acc) ^ 4 j^^syn) ^ ^ 3 s, = 0.7 m. At-*'' = 3.5 m/s. 

In the SA model with function V^7'^\g) El, 

we use 

= 23.61 m/s (85 km/h), iff "^^^ = 0.95 s'^, K^t"^ = 0.75 s^^, 
^(acc)^^-) = 0.3 + 0.4min(l,w/12) s"!, = 20 m/s, 
c - 0.007 m-i, ri^") = 1 s, = 1 m, Aw^^^ = 2.5 m/s. 



4. Diagram of Congested Traffic Patterns at On-Ramp Bottleneck in 
ATD-Model 

Numerical simulations of the ATD-model show that congested patterns (figure ^ , 
which appear on the main road upstream of the bottleneck, are qualitatively the 
same as those for the stochastic models of Ref. [SS] EHl ISH reviewed in the book [7] . 
However, dynamics of phase transitions leading to congested pattern formation and 
a diagram of these patterns in the flow-flow plane with co-ordinates are ^on and q-m 
(figure 13 (a)) exhibit some important peculiarities in comparison with the stochastic 
models P5 ,36,37. These peculiarities are associated with a deterministic character 
of the ATD-model. To understand this, firstly features of an F^S transition at the 
bottleneck in the deterministic ATD-model should be considered. 



4.I. Local Perturbation and F-~>-S Transition in Free Flow at Bottleneck 

Vehicle merging results in a abrupt local space gap reduction on the main road. This 
can lead to abrupt local vehicle deceleration. For this reason, a dynamic decrease in 
speed (figure^ (a)) and the associated increase in density in the on-ramp merging 
region appear. This local disturbance in the speed and density localized at the 
bottleneck can be considered a time-dependent dynamic perturbation in free flow. 
The dynamic nature of this perturbation (there are no random fluctuations in the 
deterministic ATD-model) is explained by dynamic rules of vehicle motion and by 
a spatial non-homogeneity localized in the on-ramp merging region within which on- 
ramp inflow and flow on the main road merge. If the flow rate qin is great enough, then 
due to dynamic merging rules of Sect. 12.61 vehicles can merge onto the main road at 
different locations within the merging region. This complex dynamic vehicle merging 
behaviour causes the associated complex dynamic spatiotemporal dependence of the 
speed and, respectively, density within the dynamic perturbation (figure 0] (a)). 

If the speed and density within the perturbation are averaged over time with an 
averaging time interval that is considerably longer than time intervals between merging 
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Figure 3. Congested patterns in the ATD-model: (a) 



Diagram of congested 
the 



patterns, (b) - Maximum capacity in free flow at the bottleneck (/maT 
discharge flow rate downstream of the congested bottleneck <?ou°'"'^'' ^^'^ 
capacity drop 5q. (c-h) — Congested patterns upstream of the bottleneck related 
to (a): (c-f) - Synchronized flow patterns (SPs) and (g, h) - general patterns 
(CPs), (c) - Widening SP (WSP). (d) - Sequence of moving SPs (MSPs). (e) - 
Localized SP (LSP). (f) ~ Alternating synchronized flow pattern (ASP), (g) - CP. 
(h) - Dissolving CP (DGP). In (c-h) the flow rates ((Jon, gin) are: (c) (350, 2140), 
(d) (60, 2367), (e) (360, 1800), (f) (240, 2026), (g) (510, 2310), (h) (360, 2310) 



vehicles/h. q^^'^°'^ ^} ; 

' ^max, lim 



2470 vehicles/h. In (b) the discharge flow rate 



changed from 2270 to 1925 vehicles/h, gout ~ 1805 vehicles/h. Tob = 30 min for 
the boundary Fi^' and 60 min for the boundary 5l^'. io = 7 min. 
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Figure 4. Perculiarities of F— >S transitions in the ATD-model: (a-f) - 
Space dependences of speed and density within local perturbations at on-ramp 
bottleneck at two different parameters of vehicle merging from on-ramp onto the 

main road Av^^^ = 8 m/s (a-c) and Av^^^ = 12 m/s (d— f). qon = Qon = 300, 

(2) 

Qon = gon = 120, and gin = 1900 vehicles/h. In (a, d) curves 1-6 are related to 
different time moments, (g) - Boundaries F^^^ of a first-order F— >S transition 

at the bottleneck associated with (a) and (d), respectively. The boundary Fg 
associated with curve 1 in (g) is taken from figured (aV In (b, c, e, f), 5-min 
averaging of traffic variables measured at virtual detectors at different locations 
are shown. Other model parameters are the same as those in figure U] 



of vehicles, then spatial distributions of the speed and density within the associated 
average perturbation (figures 01 (b) and (c)) can be considered a "deterministic" 
perturbation localized at on-ramp bottleneck. At this time scale the deterministic 
perturbation is motionless, the total flow rate (across the main road and on-ramp 
lane) within the perturbation does not depend on spatial co-ordinate. This total flow 
rate in free flow is Qsum = ftn + '/on- In contrast, the average speed and density spatially 
vary in free flow at the bottleneck. In particular, gsum = ^freipfree = v'^^^''^'' p'^^^°^\ where 
'^free Pfree minimum speed and maximum density within the deterministic 
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perturbation, respectively; the speed and density downstream of the 

perturbation, respectively (figures ^(b) and (c)); v^fj^ < v^^''^°\ p[f^^ > p(^'°®\ 

At a given gin, the greater gon, the lower the speed v^^^J^ and the greater the density 

Ph-cc within the perturbation, i.e., the greater the amplitude of the deterministic 
perturbation (figures^ (b) and (c)). This growth in the perturbation amplitude has a 
limit associated with an F^S transition that occurs spontaneously at the bottleneck 
when (Jon gradually increases. The multitude of the flow rates qin and Qon, at which 
the F^S transition occurs, determines the boundary in the pattern diagram 

(figure 13 (a)). At the boundary Fg^^ a first-order F^S transition (see Sect. I5.l|l 
occurs spontaneously during a chosen time interval Tob that is considerably longer 
than a time interval 7"^^^°^^^ (about 60 s) required for the average speed to decrease 
from the speed within a dynamic perturbation in free flow at the bottleneck to a 
synchronized flow speed (see explanations in Sect. 5.3.7 of [7])- The necessity of 
the time interval Tob is associated with a time delay Tpg for an F— ^S transition 

fB) 

found in the ATD-model: After the time delay Tpg , a time-dependent (dynamic) 
perturbation (flgure^l (a)), which can cause a short-time decrease in the speed within 

fB) 

the perturbation markedly lower than v^^J^, can occur. This perturbation occurrence 

fB) 

leads to the F— s-S transition. The boundary Fg ' is determined from the condition 

-lob ~ -^FS ■ 

In stochastic models |7], the boundary Fg^-* is also determined by the considition 

.(B) 



S 

that an F^S transition occurs at given gin and qon after a time delay T^g' during a 



chosen time interval Tob- However, in the stochastic models Tpg is a random value: 

f Bl 

In different realizations made at the same and Qon various Tpg are found. This 
stochastic model nature enables us also to calculate the probability for F^S transition 
occurrence PSIIT]. 

In contrast with the stochastic models [7], in the deterministic ATD-model there 
are no random fluctuations. Time-dependent perturbations in free flow localized at 
the bottleneck (figure^ (a)) have dynamic nature explained above. For this reason, in 

f Bl 

the ATD-model T^g' is a /jied value at given gin and Qon', consequently, the probability 
for F^S transition occurrence cannot be found. 

In addition, numerical simulations of the ATD-model show that a duration of a 

f B ) 

dynamic speed decrease within the perturbation below the speed v^^J^ is considerably 
shorter (1-3 s) than T^^ote^^^ As a result, it is found that at a given gin the time delay 

fB) fBl 

Tpg is a strong decreasing function of gon in a neightborhood of the boundary Tg : 

(B^ fBl 

Already a small increase in gon behind the boundary Fg leads to a decrease in Tpg 
down to Tj^^tJ^^n^^- Thus, we can suggest that in the ATD-model the boundary Fg^'' 
is very close to the boundary for the deterministic F^S transition (see explanation of 
the deterministic F— >S transition in Sect. 5.3.7 of Ref. [7j). 

The dynamic character of perturbations at the bottleneck, which is responsible 

f B) 

for the above mentioned physics of the boundary Fg for an F^S transition in the 
ATD-model, can clear be seen, if smaller disturbances in speed and density occur 
due to vehicle merging. Smaller disturbances can be simulated by an increase in the 
parameter Aw^^^ of vehicle merging (Sect. I2.6|l . As a result, at the same g^n and 
gon as those in figures^ (a)-(c) both time-dependent (figure^ (d)) and deterministic 
perturbations (figures^ (e) and (f)) become smaller. This leads to a shift of the 



Deterministic approach to microscopic three-phase traffic theory 



19 




Figure 5. Perculiarities of S^J transitions in the ATD-model: (a, b) — 
Spatiotemporal decay of dynamic speed waves that emerge in the merging on- 
ramp region during their upstream propagation through a widening SP (WSP) in 
space and time (a) and at different virtual detectors (b). (c, d) — Spatiotemporal 
growth of dynamic speed waves that emerge in the merging on-ramp region during 
their upstream propagation through an initial WSP in space and time (c) and at 
different virtual detectors (d). gon = 300 (a, b), 360 (c, d), and = 2310 
vehicles/h. Other model parameters are the same as those in figure 1^ 

boundary F^^^ in the diagram of congested patterns to greater gon (curve 2 in figure 01 
(g)): 

4-2. Perculiarities of S^J Transitions and Congested Patterns 

In the ATD-model, moving jam formation in synchronized flow (S— >J transition), 
which occurs at the boundary S'j^' in the congested pattern diagram (figure 13 (a)), 
exhibits also some qualitative different features in comparison with the stochastic 
models [7]- 

As in the stochastic models [7], in the ATD-model after a synchronized flow 
pattern (SP) occurs upstream of the bottleneck due to F— >S transition at the 
bottleneck, a further increase in gon leads to a subsequent decrease in the speed within 
the SP. This can cause an S-^J transition with the following general pattern (GP) 
formation. In the stochastic models, a self-growth of random model fluctuations is 
mostly responsible for the S— >J transition. In contrast, in the ATD-model there are 
no random model fluctuations. 

In the ATD-model, dynamic merging of vehicles from the on-ramp lane onto the 
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main road can cause dynamic speed and density waves that propagate upstream in 
synchronized flow of the SP (figure [SJ. It turns out that if the flow rate Qon is related 
to a point (^on, qin) between the boundaries F^^^ and (figure |2| (a)), then these 
dynamic waves decay during their upstream propagation within synchronized flow of 
the SP (figures El (a) and (b)). In contrast, at the boundary S'j^'' the waves begin to 
self-growth in their amplitude leading wide moving jam formation, i.e., one of GPs 
appears upstream of the bottleneck (figures El (c) and (d)). 

As in the KKW cellular automata (CA) model jSnHZ], in the ATD-model the 
maximum flow rate in free flow downstream of the bottleneck gmax is a decreasing 
function of ^on (figure (b)). Recall, that the flow rate qmax ^\qon) is the flow rate 
in free flow downstream of the bottleneck associated with the boundary Fg^^ . After 
a congested pattern is formed at the bottleneck, the flow rate downstream of the 
congested bottleneck called discharge flow rate gout*''"^ (figure|21(b)) is usually smaller 
than the initial flow rate gmax" The difference (Sg(gon) = qma-^ ^''('7on)-gout"'''^(9on) 
called "capacity drop" is an increasing function of Qon at the boundary Fg ' in the 
diagram of congested patterns. 

In accordance with empirical results [7|i i^ the ATD-model moving jams do not 
emerge spontaneously in free flow. This is because in all states of free flow critical 
perturbations required for an F^S transition are considerably smaller than those for 
F— >J transition. In the model, all synchronized flow states that are above the line J 
in the flow-density plane (flgure^ (c)) are metastable ones against wide moving jam 
emergence. 



5. Phase Transitions and Congested Patterns in SA-Models 

5.1. Nucleation and Metastability Effects of Pattern Formation 

As the ATD-model, the SA-niodels exhibit a flrst-order F— >S transition at the 
bottleneck, which is accompanied by nucleation and metastability effects, as well as 
by a hysteresis in SP emergence and dissolution. To illustrate these effects found for 
both the ATD- and SA-models, we restrict a consideration to the SA-model 1)35(1 - 
H41|l (flgures El and EJ . When an initial state at the bottleneck is free flow in which 
Qin is given and qon increases gradually, then, as in the ATD-model (flgure^l (a)), a 
dynamic disturbance in free flow localized at the bottleneck appears spontaneously. A 
time averaging of spatial speed and density distributions within the perturbation leads 
to the associated deterministic perturbation (flgurcs jSl (a) and (b)). Deterministic 
perturbation features are the same as those for the ATD-model (Sect.^Tj- 

The speed v\^J^ within the deterministic perturbation decreases when qon increases 
(from point 1 to 5 in flgureEl(c)). Consequently, Pf^^ increases. In the flow-density 
plane, the flow rate on the main road associated with this density increases too (from 
point 1 to 5 in flgurejSl (d)), whereas the flow rate upstream of the perturbation is 
equal to q^^, i.e., it does not change (from point 1 to 5 in flgureEl(e)). This increase 

in the deterministic perturbation amplitude with gon has a limit q^n = q'^m^'^™' '^^^ 
associated with the deterministic F— >S transition (dotted down-arrow in flgureEl(c)). 
However, non- homogeneous free flow dynamics fSect. ITT|l . which is caused by vehicle 
merging, results in an F— >S transition at a smaller qon (point 5 in flgureEl(c)) related 
to a point on the boundary F^'^ in the diagram of congested patterns (flgurejSl (a)). 
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Figure 6. Nucleation and mctastability of free flow and synchronized flow: (a, 
b) — Spatial dependence of the speed (a) and density (b) within the deterministic 
perturbation in free flow at the bottleneck for two flow rates qon = 150 (curve 
I) and 240 (curve II) vehicles/h. (c) - Double Z-characteristic in the flow-rate- 
speed plane for the sequence of F— >S— >J transitions, (d, e) — Hysteresis effects 
in the flow— density plane due to F— >S and reverse S— >F transitions within the 
deterministic perturbation (d) and on the main road upstream of the bottleneck 
(e). Dotted curves in (d, e) are related to steady state model solutions, gin = 
2252 vehicles/h. 
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Figure 7. Congested pattern excitation in metastable free flow: (a) - Short-time 
(2 min) perturbations in tlic flow rate qon used for pattern excitation, (b— d) - 
WSP (b), GP (c), and MSP (d) induced by tiie related perturbations (curve 1, 2, 
and 3 in (a), respectively), gin = 2252 vehicles/h. qon = 220 vehicles/h in (b, c) 
and Qon = 140 vehicles/h in (d). Amplitudes of perturbations in (a) are: 270 for 
(b, d) and 500 vehicles/h for (c). 



The speed decreases and density increases abruptly within the initial perturbation 
(arrows F— >S from point 5 to 5' in figures|ni(c)-(e)) and a congested pattern emerges 
at the bottleneck. In the example, a widening SP (WSP) occurs upstream of the 
bottleneck due to the F^S transition (figure El (c)). 

If now gon decreases, the speed within the WSP increases (from point 5' to 3' 
in figures (c)-(e)). This synchronized flow speed increase has a limit: The speed 
increases and density decreases abruptly within the synchronized flow (arrows S— >F 
from the point 3' to 3 in figures El (c)-(e)) and free flow returns at the bottleneck. 

Note that qi^ in flgure is chosen to be greater than the threshold flow rate gth 
for moving SP (MSP) existence. As a result, the initial motionless downstream front 
of synchronized flow at the bottleneck begins to move away upstream. Consequently, 
an MSP emerges (figure (HI (d)) (range of gon within which MSPs occur is shown by 
a dashed part of the synchronized flow states Wsyn in figure El (c)). At greater ^on 
on the dashed part of the synchronized fiow states Vsyl in figure El (c) this free flow 
at the bottleneck can persist for a short time only: A new F— >S transition occurs 
spontaneously and a new MSP emerges at the bottleneck, and so on. Due to this 
effect, a sequence of MSPs appears. 

At qin > qth an MSP can also be induced by application a short-time local 
perturbation in free flow. The speed within this external perturbation should be lower 
than the critical speed v^^ pg associated with the critical branch on the Z-characteristic 
for the F— »S and reverse S— »F transitions at the bottleneck. As in the stochastic 
models |7], this Z-characteristic consists of the states for free flow associated with 
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Figure 8. Diagram of congested patterns at the on-ramp bottleneck in SA-model 
JSg-iHJ (a), the maximum capacity in free flow Qmax^ ^ the discharge flow rate 
at the on-ramp and the capacity drop Sq (b), and congested patterns 

(c-h) related to (a): (c-e) - SPs and (f-h) - GPs. (c) - WSP. (d) - MSP. (e) 

- LSP. (f) - GP arising from WSP at smaller Qon- (g) - GP at gin > gout- (h) 

- GP at gin < gout. In (c-h) the flow rates {gon,gin) are: (c) (200, 2353), (d) 
(140, 2378), (e) (940, 1520), (f) (250, 2353), (g) (400, 2353), and (h) (900, 1770) 
vehicles/h. n^^^'^'^ ^) 2475 vehicles/h. 

' ^max, lim ' 
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the deterministic perturbation at the bottleneck v^^J^, the critical branch v^^ pg, and 

synchronized flow states Vsyn (figure |d1 (c)). In accordance with this Z-characteristic, 
we get the associated hysteresis effects on the fundamental diagram (arrows F^S and 
S^F in figures ini (d) and (e)). 

If in contrast gon increases, the speed within the WSP decreases (figure El (c)). 
This speed decrease has a limit associated with the flow rate gon — qin' at which 
an S— >J transition must occur (dotted down-arrow S— in figure El (c)). However, 
because there are speed and density waves of a finite amplitude in synchronized flow, 
an S— s-J transition occurs already for gon < 9ori'' ^"^^ (point 6 and solid down-arrow 
in flgure El (c)). As a result, an GP emerges. This is because in the SA-models, 
synchronized flow steady states with the speed v > Vcr , which are above the line J 
in the flow-density plane, are metastable ones against wide moving jam emergence. 
This metastability can be seen from another Z-characteristic in the speed-flow plane 
associated with an S— >J transition in synchronized flow. The Z-characteristic consists 
of the states for synchronized flow Vsyn, the critical branch for critical perturbations 
in synchronized flow t^^f^sj' ^^'-^ ^^^^ w = for wide moving jams (flgureEl (c)). 

From the resulting double Z-characteristic (flgure El (c)), it can concluded that 
in a metastable free flow at the bottleneck (left of the boundary Fg in the diagram 
in flgure IHI (a)) depending on amplitude of a time-limited perturbation caused, for 
example, by an increase in qon (curves 1 and 2 in flgure|Z|(a)), either an WSP (flgure|Z| 
(b)) or an GP (flgure [7| (c)) can be induced. At smaller Qon (curve 3 in flgure [7| (a)), 
an MSP (figure |7| (d)) can be excited in free flow. All results presented in figures El 
and|7|for the SA-model remain qualitatively equal for the ATD-model. 



5.2. Comparison of Congested Patterns in ATD- and SA-Models 

The SA-model (|36f) - (|41|l (figure EJ exhibits the following shortcoming in comparison 
with the ATD-model (figure Ell: 

(i) If the flow rate qon is within a flow rate range gon^ '"' < Qon < gon^^' , then no 
SP can be formed at the boundary Fg ' in the diagram (flgure |8| (a)): The sequence of 
F— *S— >J transitions occurs spontaneously at this boundary, leading to GP emergence. 
For this reason, the related part of the boundary at which GPs emerge spontaneously 
in free flow at the bottleneck is labelled & 'S'j • 

(ii) If the flow rate gin at this boundary decreases, another characteristic flow rate 

fLSPl . . fLSPl 

gin = gin associated with the flow rate gon = gon at this boundary is reached: 
At gin < gin moving jams do not emerge in synchronized flow upstream of the 
bottleneck. As a result, at gin < 9in and right of the boundary Fg^^ only an LSP 
remains at the bottleneck. Within this LSP the speed is very low. This LSP has a 
qualitative different nature in comparison with an LSP of higher synchronized flow 
speed in the ATD-model that occurs at considerably greater gin (figure EJ. 

In the SA-model the branch for average synchronized flow states 

Vav^'^^ has a part with a positive slope (flgureEl (d)). Then LSPs of higher speeds 
appear in the diagram of congested patterns (flgure EJ. However, these LSPs are not 
related to LSPs observed in empirical observations. To explain this, note that these 
model LSPs are very narrow ones (flgureEl (e)). They are localized within the merging 
region of the on-ramp and consist of two narrow fronts only fflgure lTUl fa')'): There is no 
region of synchronized flow between the fronts within these LSPs. This is regardless 
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of the flow rates q-m and gon- Conflictingly, in empirical observations rather than 
such narrow LSPs, an extended region of synchronized flow is usually observed within 
an empirical LSP. The LSP width (in the longitudinal direction) changes over time 
considerably. These empirical features of LSPs shown by the ATD-model (figure E| 
(b)) are not found in the SA-model. 

(iii) In the ATD-model (figuresEH (e) and (f)) as in empirical observations, both 
free and synchronized flows can be formed between wide moving jams within an GP. 
In contrast, in the SA-models only free flow can be formed between wide moving jams 
within the GP (figuresEl(c) and (d)). The reason for this is as follows. The average 
branch for synchronized flow lies for speeds v > Wcr'''' above the line J (figure|21 (g))- 
Flow states in the jam outflow should be related to points on the line J. Thus, there 
are no synchronized flow states between the jams. This explains why only free flow 
can be formed between the jams in the SA-models.* 

The mentioned shortcoming of the SA-modcl result from the averaging of a 2D- 
region of steady states for synchronized flow in the flow-density plane of the ATD- 
model (figure n (a)) to the branch for average synchronized flow states (curve S in 
figure EJ. 

At chosen SA-model parameters the condition 

gout > g^P'"^'^) (48) 

is satisfied, where g(P™'^'') is the fiow rate within the pinch region of an GP in which 
narrow moving jams emerge. Under the condition (|48|l . no DGPs appear in the SA- 
models (figures|Hland|ni). At other parameters of the SA-models, an opposite condition 

gout < g^P'"^''^ (49) 

can be satisfied. Then DGPs appear in the SA-models. 

The maximum flow rate in free flow downstream of the bottleneck gma'x ^^(gon), 
the discharge flow rate gout**''^\ ^.nd the "capacity drop" Sq can sometimes exhibit 
different features as those in the ATD-model (figure 13 (b)) when the flow rate gon 
changes (figures |S1 (b) and|ni(b)). Particularly, in contrast with the ATD-model, in 
the SA-model Ip!5|l - (|1T|| gmax does not depend on gon, whereas in the SA-model 
lp!^ - lP?1| . (|^ gmax depends on gon but at considerably greater gon than for the 
ATD-model. 

Simulations show that the SA models presented in |Appendix B| show qualitatively 
the same features of the phase transitions and spatiotemporal congested patterns as 
those in the SA-model if^ - ljlTT) . 



5.3. Comparison with Stochastic SA-Models 

It is interesting to compare the deterministic SA-models with possible stochastic SA- 
models. Such models can be derived from the stochastic model of Ref. |35[ I37| . if 

* The only exclusion is the SA-model with the average branch for synchronized flow I42i . if 
parameters for the curve S and/or the line J, i.e., for wide moving jam propagation are chosen 
different as those shown in figures |21 (c) , (d), and (h): These different parameters should lead to an 
intersection of the line J with the average branch for synchronized flow with a positive slope in the 
flow-density plane. However, in this specific case only one state of the synchronized flow, which is 
associated with the point of the latter intersection, is possible. This model effect is not agreed with 
empirical results, in which the flow rate and speed between wide moving jams within GPs can change 
over time considerably |7|. 
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Figure 9. Diagram of congested patterns at the on-ramp bottleneck in the 
SA-model 1351 — EHl . I42i (a), the maximum capacity in free flow gmax ^\ the 
discharge flow rate at the on-ramp 90^°**'°' and the capacity drop 5q (b), and 
congested patterns (c-f) related to (a): (c-e) - SP and (f) - GP. (c) - WSP. (d) - 
Subsequence of MSPs. (e) - LSP. (f) - GP arising from WSP. In (c-f), the flow 
rates ((jon,gin) are: (c) (200, 2400), (d) (90, 2450), (e) (360, 2115), (f) (300, 2400) 
vehicles/h. q^^^ ~ 2475 vehicles/h. In (b), is changed from 2450 to 

2200 vehicles/h. 'gout ^ 1880 vehicles/h. 



2D region of synchronized flow steady states is averaged to synchronized flow states 
related to a ID region in the flow-density plane. 

A stochastic SA-model can easily be derived from the stochastic model of Ref. I37| 
based on the physics and ideas for the SA-model approach discussed in Sect. 13 To 
reach this goal, in the part of the stochastic model of Ref. |37] 

Vn+l = max(0, min(wfroo, fc,n, Ws,n)), (50) 
Xn+1 = Xn + Vn+lT (51) 

for a desired speed in synchronized flow Wc,nj rather than the formula (3) of Ref. liJTj 
leading to a 2D region of synchronized flow steady states in the flow-density plane, 
the following equations associated with the physics of the SA-models of Sect. Oare 
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synchronized flows 




Figure 10. Comparison of LSPs (a, b) and CPs (c-f) in the SA-model 
I42i (a, c, d) and in the ATD-modol (b, o, f). LSPs and their parameters are the 
same as those in figures |5] (e) and|3]{e), respectively. In (e, f), Ai;^^' = 7.5 m/s. 
(9on, gin) are: (c, d) (380, 2397), (e, f) (660, 2222) vehicles/h. 



used: 

Vc,n = Vn + max(-6„T, niin(a„T, A„)), (52) 
_ . K{vtn - t;„) at g > g^^^^^ , 

In (|50|) - (|53|l . and Xn are the speed and space co-ordinate of a vehicle; the index 
n corresponds to the discrete time t — nr, n = 0,1,2,..; r is the time step; Wfi-eo 
is the maximum speed in free flow, which is a constant; Vs^n is the save speed of 
Ref. On > is acceleration, 6n > is deceleration, which are taken as the same 
stochastic functions used in the model of Ref. jSZ]; the space gap gn = xi^n — Xn — d; 
the average speed in synchronized flow steady states 14v'*"^'' is given by the formula 
l)41|l at gma™'' = 0. Of course, other formulations for the average synchronized flow 
steady states Vav^'^\ for example used in the deterministic S A- models (figures [21 (b), 
(d), and (f)) can also be applied. 

In general, descriptions of random vehicle acceleration and deceleration are the 
same as those in the stochastic model of Ref. [SSJE?]: At the first step, the preliminary 
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speed Vn+i is set to Vn+i — ^n+i where the speed is calculated from the equations 
l(Sn|l - lf5^ . At the second step, a noise component is added to the calculated speed 
■Un+i. Then the final speed is found from the condition |37| 



where Omax is the maximum acceleration. 

However, in contrast with the stochastic model of Ref. |37], the noise component 
in 154|l is chosen to be different from zero only if the vehicle decelerates, specifically 



where r = rand(0, 1), 9{z) = at z < and 9{z) = 1 at z > 0, (5 ^ TOmax, ^max, Ph 
are constants. 

Simulations show that the stochastic SA-model H50 |l - H55|l exhibits qualitatively 
similar spatiotemporal congested patterns at the on-ramp bottleneck (figure 111(1 as 
those in the associated deterministic S A- models (figure IS)). However, there are 
qualitative differences in the dynamics of first-order F— +S and S— >J transitions leading 
to pattern formation explained in Sect. 15. l1 In the stochastic SA-model, random model 
fluctuations are important for phase transition nucleation, whereas in the deterministic 
SA-models the F— >S and S—^J transitions are nucleated by dynamic perturbations 
emerging within the on-ramp merging region. 

Note that under the chosen model parameters in the stochastic SA-model (|50fl - 
H55|l the condition H49|l can be satisfied at smaller Qon- As a result, there is a region 
in the diagram of congested patterns in which DGPs occur (region labelled DGP in 
figure im fa)). After the wide moving jam of the DGP is upstream of the bottleneck 
as well as in the ATD-model, an LSP remains at the bottleneck ffigure ITTl (h)). In 
contrast with the ATD-model, this LSP exists for a finite time interval only and free 
flow returns at the bottleneck. 

However, in a small neighbourhood of the boundary labelled D in the diagram, 
which separates DGPs and GPs, there is a peculiarity in pattern formation under the 
condition H49|l : If the flow rate gon increases in the neighbourhood of the boundary D 
in the diagram, then the lifetime of an LSP, which occurs within an DGP increases 
and it tends towards infinity at the boundary D (figure El (a)). This quasi-steady 
LSP is explained by a very long interval between wide moving jam emergence in the 
synchronized flow at the bottleneck (flgure^l (b)). This interval tends towards the 
infinity at the boundary D. 

6. Discussion 

6.1. Comparison of ATD- and SA-models with OV-models and other Deterministic 
Models 

The first term in the formula for vehicle acceleration 5^^''°°^ (|12|l . ^(^^^''""^(g) — v), 
describes the dynamics of the speed v in the vicinity of the optimal speed V^^'"^''\g) 
in free flow. At a time scale that is considerably greater than the time delay r, this 
dynamic behaviour is the same as those in different OV-models E El 121] , which can 
be written as follows 



Vn+l = max(0, min(wfroc, Vn+l + Cn, Vn + ama^T, Vs^n)), 



(54) 




(55) 



- = A{g, v){V{g)~v). 



(56) 
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Figure 11. Diagram of congested patterns at the on-ramp bottleneck in the 
stochastic SA-model l5Ui - l55l (a), the maximum capacity in free flow gmaT ^\ 
the discharge flow rate at the on-ramp and the capacity drop 5g (b), 

and congested patterns (c— h) related to (a): (c— f) — SPs and (g, h) — GPs. 
(c) - WSP. (d) - MSP. (e) - LSP. (f) - ASP. (g) - GP. (h) - DGP. In (c-h) 
the flow rates (gon,gin) are: (c) (200, 2250), (d) (60, 2250), (e) (720, 1470), 
(f) (120, 2235), (g) (500, 2220), and (h) (352, 2235) vehicles/h. q^^^ fjj^ « 

2300 vehicles/h. Model parameters: A(f"'°<=) (3,1) = 0.5 min (l, {g^ - )/20) , 



g^;^") = 36 m, = 0.1, T^^^^ = 1.45 s, g^ax^ = 0, Pb = 0.02. pi =0.55, 

P2{vn) = 0.5 + 0.48e(Dn - 15), amax = bmax = 0.5 m/s^, K = 1. The other 
parameters are the same as those in |33. The on-ramp model from 1371 1381 is 

used, in which the synchronization gaps G„ = G„ = g ■ ■ 
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Figure 12. Transformation of DGP in figure in]{h) to GP with a very long 
time interval between wide moving jam emergence (b) through an DGP with a 
quasi-steady LSP (a) by a small increase in the flow rate gon in a neighbourhood 
of the boundary D in the diagram of congested patterns in figurc lTTI fa'l. The flow 
rates (ijon.gin) are: (a) (357, 2235) and (b) (362, 2235) vehicles/h. 



However, in (|56|l the vehicle acceleration A{g, v){V{g) — v) is valid for the whole 
possible space gap range [21^13 

.9 > 0. (57) 

In contrast with the OV-models, in the ATD-model this vehicle acceleration is applied 
for large space gaps |SJ associated with free flow only. 

The crucial difference of the ATD-model with the OV-models and all other 
deterministic microscopic traffic flow models (see references in the reviews 
is that the vehicle acceleration behaviour qualitatively changes when the vehicle is 
within the synchronization gap, i.e., if the condition I0J), which is opposite to the 
condition Q, is satisfied. 

The condition is associated with the synchronized flow phase in which there 
is no optimal speed in the ATD-model. This conclusion follows from (I13II and its 
analysis made in Sect. 12.81 in which it has been shown that for a given steady space 
gap in synchronized flow there are an inflnity of steady vehicle speeds within a finite 
speed range (flgure^ (a-))- 

The concept of safe speed Vs{g, Vi) for vehicle collision prevention used in the 
ATD-model is qualitatively different from the concept of optimal speed that is the 
basis of the deterministic approaches (01 and Q : The optimal speed is a desired one 
(this explains the term "optimal" speed) for a driver to be reached (the driver moves 
comfortable with the optimal speed during a long time) , whereas the safe speed is not 
an optimal one but a limiting speed that is still permitted (the driver should not move 
with this speed during a long time because this is strain for the driver and, therefore, 
non-comfortable). The qualitative difference of these two concepts is mathematically 
reflected in the dynamic model behaviour. In the ATD-model, when the vehicle speed 
is higher than the safe speed and safe deceleration is applied, then a driver time delay 
is equal to a small driver reaction time: t — Ts H21|I . In all other driving situations, 
which are not associated with safe speed, driver time delays are different from Tg. This 
is because these driver time delays are associated mostly with qualitatively different 
expected events occurring within different traffic phases (Sect. IT^ . As a result, in 
the ATD-model driver deceleration to the safe speed occurs considerably quicker, then 
in other driving situations. In contrast, in accordance with the concept of optimal 
speed, in OV-models there is a driver time delay in deceleration that characterizes 
speed relaxation to the optimal speed P H E3 EHl EHl EOI ■ 
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The crucial differences between the SA-models and ah other traffic flow models in 
which steady states covering a one-dimensional region(s) in the flow-density plane (see 
references in the reviews are as follows. In contrast with the models of 

Ref. U El El El El J in the SA-models at each density of free flow states the critical 
amplitude of a local perturbation required for an F— >S transition is considerably 
smaller than the critical amplitude of a local perturbation required for an F^J 
transition. 

In the SA-models, there are two ranges of model steady states separated one from 
another by a model discontinuity in vehicle space gap or in speed (figures |21(a)-(e)) or 
else due to instability of model steady states against infinitesimal non-homogeneous 
fluctuations (figure El (f)) This simulates the hypothesis of three-phase traffic theory 
about a competition between over-acceleration and speed adaptation effect that is 
responsible for F— >S and S^F transitions: The first range of steady states simulates 
free flow, whereas the second simulates synchronized flow. To simulate S— >J transitions 
within synchronized flow, steady states associated with synchronized flow of higher 
speeds are metastable with respect to moving jam emergence, i.e., moving jams emerge 
in these synchronized flow states only if large enough amplitude local perturbations 
appear; synchronized flow states of lower speeds are unstable with respect to moving 
jam emergence. These requirements to the SA-models lead to F— >S^J transitions 
that are responsible for moving jam emergence found in empirical data 0. 



6.2. Critical Discussion of Theories and Models based on the Fundamental Diagram 
Approach 

In the OV model H56|l . as in other deterministic (and stochastic) traffic flow models 
in the context of the fundamental diagram approach reviewed in ^ El El El EI > which 
claim to show spontaneous jam emergence, there is a range of the density on the 
fundamental diagram in which steady states on this diagram are unstable against 
infinitesimal perturbations. U This instability leads to wide moving jam emergence 
in these models both on homogeneous road and at a bottleneck. We denote the 
minimum density of this density range, in which infinitesimal fluctuations grow, by 
pii'' (figures US (a) and (b)). 

There are two possibilities for the arrangement of the point of this instability 
(pcr\ Qc?) on the fundamental diagram in the OV model and other models in the 
context of the fundamental diagram approach: 

(i) The point {pii\ Icr^) lies left of the maximum point of the fundamental 
diagram (po, Qo), i-e., on the branch of the diagram with a positive slope f figure [T^ 

(a) ): 

P = pi'? < Po (58) 

(ii) The point {pii\ qii"^) lies right of the maximum point of the fundamental 
diagram (po, Qo), i-e., on the branch of the diagram with a negative slope f figure [T51 

(b) ): 

P = pi? > PO. (59) 

tt It should be noted that some of the models based on the fundamental diagram approach are not 
valid far from equilibrium. It is not simply a matter of a phase transition type that a model exhibits in 
steady conditions, but of the difficulty of closing equations, which should work in unsteady conditions 
by relations valid only in steady uniform conditions. 
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Figure 13. Traffic patterns in OV-models at an on-ramp bottleneck: (a, b) - 
Fundamental diagrams of OV-models when the condition 1581 (a) and I59i (b) 
are satisfied, respectively, (c, d) — Widening patterns of dense flow upstream 
of the bottleneck at two flow rates (jon =160 (c) and 400 (d) vehicles/h. (e) - 
Formation of wide moving jams within the dense flow upstream of the bottleneck 
at qon =700 vehicles/h. In (c-e) the flow rate =2676 vehicles/h. Figures (c— e) 

are related to the OV- model (b) at V{g) = Vb ( tanh((g — go) /gi) + tanh(go /gi)) 
at Vb = 14 m/s, go = 17 m, gi = 7 m, the sensitivity A{g, v) = A{v) in I56i is 
A(v) = 5 at ?; > 12 m/s and A{v) = 0.9 at v <12 m/s. 



Note that in both cases (i) and (ii) all states on the fundamental diagram, in 
which the density satisfies the condition 

Pn.in<P<pii\ (60) 

where pmin is the density in the wide moving jam outflow associated with the flow 
rate gout, are metastable states with respect to moving jam emergence The case 
(i) has intensively been considered in the literature and criticized in Sect. 

3.3.2 of the book [7j. 

In the case (ii) (figureE|(b)), the flow rate in free flow downstream of an on-ramp 
bottleneck gsum cannot exceed the maximum flow rate on the fundamental diagram 
qo . If (7in is a given large enough value and the flow rate gon begins to increase, then a 
localized perturbation as that in the ATD- and SA-models (Sects. ETI and l5!T|) appears 
at the bottleneck (figures (a) and (b)). The minimum speed within the time- 
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Figure 14. Pattern features in tlie OV-model in figure 1131 (b) at an on-ramp 
bottleneck: (a, b) - Spatial dependences of the averaged speed (a) and density 
(b) on the main road within deterministic perturbations localized at the bottleneck 
at qou = 80 (curve I) and 117 (curve II) vehicles/h. (c, d) - On-ramp flow rate 
dependencies of the average speed (c), flow rate and density (d) on the main 
road at locations of the minimum of the average speed, (e) - On-ramp flow rate 
dependences of the flow rate and density on the main road at the location 200 m 
upstream of the begin of the on-ramp merging region, (f, g) - On-ramp flow rate 
dependences of the flow rate (f) and density (in the flow— density plane) (g) on the 
main road in free flow downstream of the bottleneck, (h) - Wave of dense flow 
that starts to propagate upstream with the velocity ^ — 0.7 km/h under the 
condition I62i at a very small value Aq = 3 vehicles/h; ti = 15, t2 = 25, ts = 55, 
t4 = 90 min; Qon — 120 vehicles/h; the deterministic perturbation (curve II) is the 
same as those in (a). Dotted curves on (d, e, g) show steady model states in the 
flow-density plane, = 2676 vehiclcs/h, qq = 2795 vehicles/h. 5-min averaged 
data. 
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averaged (deterministic) perturbation on the main road decreases when increases 
(points 1-3 in figure H"!] (c)). In the OV-model, when ^on increases, the location of 
the minimum speed within the deterministic perturbation on the main road exhibits 
firstly a shght shift downstream and then upstream within the merging region of the 
on-ramp (figure ITU (a, b)). For this reason, the flow rate on the main road at the 
location of the minimum speed within the deterministic perturbation on the main 
road firstly slightly increases and then decreases (points 1-3 in figure El (d)), whereas 
the flow rate on the main road upstream of the on-ramp merging region is equal to q-m 
(points 1-3 in figure ITU (e) ). When qon increases beginning from zero, the flow rate 
Qsum downstream of the bottleneck increases beginning from (?sum = 9in (points 1-3 in 
figures [TH ff ) and (g)). 

At a given large enough flow rate gin this growth of the local perturbation in free 
fiow at the bottleneck with ^on has a limit. This limit is reached when the flow rate 
Qon reaches some critical value gon = gon^ at which the flow rate gsum is equal to the 
maximum flow rate on the fundamental diagram: 

gsum = gin + g^? = go- (61) 

When the flow rate gon increases further, i.e., 

Ag gin + gon - go > 0, (62) 

then the upstream front of the initial perturbation, which is motionless at the condition 
gsum = gin + gon < go (curve II in figure 1711 (h)), begins to move upstream of the 
bottleneck, i.e., a wave of lower speed and greater density propagating upstream 
appears (spatial speed distributions related to the times ti-t^ in figure ITU fh)). As a 
result, a dense flow associated with the branch of the diagram with a negative slope 
occurs upstream of the bottleneck (figures [TSI (c) and (d) and points 4-6 in figures ITU 
(c)-(e)). At the critical point (|61|l . the derivative of the minimum average speed on 
the main road on the fiow rate gon is discontinuous, whereas this speed is a continuous 
decreasing function of gon (figureE|(c)). The greater the flow rate gon, specifically, the 
greater Ag (|62f) , the greater absolute velocity of the wave of dense flow propagation 
\vd\ (figures El (c) and (d)). In addition, the flow rate downstream of the bottleneck, 
which is equal to go under the condition (|61|l . remains approximately to be equal to 
go, when gon increases (points 4-6 in figures IT^ ({) and (g)). 

It must be noted that the above mentioned behaviour of the upstream front of 
the perturbation at the bottleneck in the OV model (|56l) , (|59|) is qualitatively different 
from those for the upstream front of the perturbation at the bottleneck in the ATD- 
and SA-models. In the latter case, when the flow rate gon reaches the critical value 
for an F^S transition, a wave of synchronized flow occurs abruptly and propagates 
upstream with a finite velocity. This is associated with a first-order F— i-S transition. 
In contrast, in the OV-model there is no discontinuous change in the velocity Vd when 
due to an increase in gon the condition (|62|l is satisfied: |ud| increases continuously 
beginning from zero, when gon first reaches and then exceeds the critical fiow rate gon"* 
associated with the condition Specifically, we find that if Ag — * 0, then |ud|^ 0. 
Thus, in the OV-model there is no first-order phase transition from free flow to dense 
flow. 

The widening dense flow upstream of the bottleneck (figuresEl(c) and (d) and 1 141 
(h)) can exist only, when the density pd in the dense fiow satisfies the condition 

Po<Pd<piV. (63) 
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This is because at the density p^i = Per the dense flow loses its stability against wide 
moving jam emergence (point 7 and dotted down-arrow in figure [TH (c)). However, 
dynamic waves that emerge due to vehicle merging at the bottleneck propagate 
through the dense flow. For this reason, in numerical simulations this moving jam 
emergence occurs already at the density < Per'' (point 6 and solid down-arrow in 
figure [T3|(c)). 

The congested patterns in figures El (d) and (c) at the first glance resemble 
a widening SP and an GP, respectively. Indeed, in both cases a dense flow occurs 
upstream of the bottleneck whose downstream front is fixed at the bottleneck. Thus, 
this dense flow should satisfy the macroscopic spatiotemporal objective criteria for 
the synchronized flow phase (Sect. This conclusion is, however, incorrect. To 
explain this, note that in empirical observations application of the objective criteria, 
which define the traffic phases in congested traffic, leads to clear distinction of the 
synchronized fiow phase. This synchronized fiow exhibits the following fundamental 
empirical feature: An F— >S transition leading to synchronized flow emergence is a 
first-order phase transition. In contrast, in a traffic flow model an application of 
the objective criteria does not guarantee that dense flow occurrence in free flow is 
associated with a first-order phase transition, which is one of the requirements for the 
synchronized fiow phase. 

This conclusion concerns the OV model (|56|l . H59|l (figure El (b)) as well as other 
models in the context of the fundamental diagram approach under condition H59|l . 
Whereas for the SA-model there is a Z-shaped speed-flow characteristic associated 
with a flrst-order F— >S transition in free flow at the bottleneck (flgures (c)-(e)), 
for the OV model the on-ramp flow rate dependence of the speed at the bottleneck 
is a monotonous decreasing function fflgure ITll (c)): There is no flrst-order phase 
transition, when a dense flow related to the fundamental diagram with a negative 
slope is formed upstream of the bottleneck. Thus, the dense traffic flow in the case of 
the OV model and other models in the context of the fundamental diagram approach 
under condition l|59|) does not exhibit the important empirical feature of synchronized 
flow and, therefore, the dense flow is not associated with the synchronized flow phase. 

There are also traffic flow models in the context of the fundamental diagram 
approach, in which there is no instability of steady model states on the fundamental 
diagram regardless of the vehicle density. Examples of this model class are as follows: 
(i) An OV model (|56|l in which the sensitivity A{g, v) is great enough regardless of 
V and g. (ii) The Nagel-Schreckenberg cellular automata model in the deterministic 
model limit, i.e., when probability of model fluctuations in this model is equal zero 
{p = 0) ■ (iii) The Lighthill-Whitham- Richards model 47 and the associated cell- 
transmission models 0Hj- In this model class, traffic patterns at a freeway bottleneck 
are qualitatively similar as those found in the OV model (|56|) . (|59|l at the density 
considerably smaller than the critical density pcr^ (the patterns associated with points 
1-5 in figuresIT^ (c')-(o')'). These common model features are as follows: 1) the local 
perturbation at the bottlenecks at Aq < ffigureslTH fa) and (6)); 2) widening dense 
flow upstream of the bottleneck at Aq > fflgures lTST c) and (d) andll4l (h)): 3) there 
is no discontinuous change in speed (no speed breakdown) at the bottleneck when 
widening dense flow occurs; 4) with an increase in traffic demand at Aq > 0, the 
upstream front velocity of widening dense flow increases continuously beginning from 
zero. Thus, in this model class, there is no flrst-order F— i'S transition observed during 
the onset of congestion at the bottleneck, i.e., this dense flow has no relation to real 
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freeway traffic. 
6.3. Conclusions 

(i) Two different deterministic microscopic traffic flow model classes in the context 
of three-phase traffic theory, the ATD- and SA-models, have been introduced in the 
article. 

(ii) The ATD- and SA-models reproduce important empirical spatiotemporal 
features of phase transitions in traffic flow and congested traffic patterns. 

(iii) In contrast with all other known deterministic microscopic traffic flow models, 
in the ATD- and SA-models vehicles moving in free flow and vehicles moving in 
synchronized flow exhibit qualitatively different dynamic behaviour. This is a result 
of the introduction of two separated regions of steady state model solutions for free 
flow and synchronized flow in the ATD- and SA-models as well as different dynamic 
rules of vehicle motion in free flow and synchronized flow implemented in the models. 

(iv) As in empirical observations, there is a first-order phase transition in the 
ATD- and SA-models from free flow to synchronized flow that explained the onset of 
congestion at bottlenecks in these models. 

(v) The nature of the onset of congestion as a first-order F— >S transition in free 
flow at the bottleneck, which the ATD- and SA-models show, is also associated with 
metastability of free flow at the bottleneck against external short-time disturbances in 
this flow in a neighbourhood of the bottleneck. As a result, there is multiple congested 
pattern emergence in an initial free flow at the bottleneck in the ATD- and SA-models: 
Depending on an amplitude (or duration) of an external disturbance, one of the SPs 
or else an GP can be induced in free flow at the bottleneck at the same chosen model 
parameters. 

(vi) In accordance with empirical results, in the ATD- and SA-models moving 
jams can emerge spontaneously in synchronized flow only, i.e., as a result of F— i'S^J 
transitions. 

(vii) In addition to the above common behaviour of the ATD- and SA-models, 
these models exhibit also some qualitatively different features. This is because in the 
ATD-model synchronized flow model steady states are related to a 2D-region in the 
flow-density plane, whereas synchronized flow model steady states in the SA-models 
belong to an ID-region (a curve) in the flow-density plane. In particular, the following 
differences of model features have been found: 

(1) The ATD-model can show all types of spatiotemporal congested patterns at 
an on-ramp bottleneck observed in empirical observations. 

(2) In contrast, SA-models cannot show LSPs associated with empirical results 
as well as some of empirical features of synchronized flow between wide moving jams 
within GPs. 

(viii) Models in the context of the fundamental diagram approach reviewed 
in m^lEl^lSlE] cannot explain the onset of congestion in free flow, which in empirical 
observations is associated with a first-order F^S transition. Depending on the model 
type and model parameters, in these models either wide moving jam emergence is 
responsible for the onset of congestion at an on-ramp bottleneck rather than an 
empirically observed F^S transition or a widening dense traffic flow occurs upstream 
of the bottleneck when the density in free flow at the bottleneck exceeds the density 
associated with the maximum point on the fundamental diagram. In the latter case, 
in contrast with empirical observations there is no first-order phase transition from an 
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initial free flow to this dense flow at the bottleneck: The dense flow results from non- 
homogeneity of a freeway in a neighbourhood of the bottleneck. Thus, these models 
cannot show a first-order F^S transition observed during the onset of congestion 
at the bottleneck in real freeway traffic, i.e., this dense flow has no relation to real 
freeway traffic. Indeed, the first-order F— >S transition is a fundamental empirical 
feature of the onset of congestion in free flow with the subsequent synchronized flow 
phase emergence at the bottleneck. 



Appendix A. 

To derive formula 1)181) |45| . let us consider a solution of (|16|l when it is an equality: 

Us q, Vf) — , A.l) 

From H16() . IjA.ip . it can be seen that if g — v^Ts, then the safe speed Vg = vg; if 
in contrast g < viT^, then the speed Vg < v^. In particular, this ensures collision 
less vehicle motion. To simplify the formula (|A.1|) . let us replace the space gap g in 
denominator of (|A.ip by the value VfTg. This reduces the safe speed Vg at g < veTs, 
therefore, the safety condition H16|l remains to be valid. Then from formula (|A.1|I . we 
get 

g + vj/{2bs) 

Vsig, vi) = — — — — . A.2) 

Ts + Vi/[2bs) 

To provide more comfortable vehicle deceleration, an anticipated gap g*^^^ = g + {ve — 
v)To is used in formula (|A.2|I rather than the gap g. As a result, (|A.2|I takes the form 

g + {v,-v)To + vj/i2h) 
Vsig, ve) = — — — — — . A.3) 

Ts + vi/[2bs) 

Substituting I) A. 3(1 into we find formula (^HJ with coefficiens (EHJ, (EOJ- 

Note that we have also tested another formulation for the safe speed in the ATD- 
model when the speed Vs{g, Vf) in l|15() is given by formula (|A.ip . Simulations of the 
ATD-model show that both formulations ((18|l and (|A.ip ensure collision less vehicle 
motion at an appropriate choice of model parameters and lead to qualitatively the 
same features of phase transitions and congested patterns. 



Appendix B. 

In this Appendix, two further variants of the SA-models are presented. In the first of 
these variants, the formula 136(1 reads as follows 

at 5 ^ ffmin ' 

^ig'i^'' <g<g^^^\ (B.l) 

^ r\ ^ ^ (jam) 
at < 5 < ffmax , 

where a^^''""), a^'^"), a^^'") are given by In this SA-model, steady states 

of free flow (the curve F in figure |3(e)) correspond to the condition ((43(1 . averaged 
steady states of synchronized flow are related to a line S given by the condition 

'/-(i-p/p|ir^)/7^ir^ ^(^<p<(^. (B.2) 



dv 
'di 
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steady states for a wide moving jam are associated with the condition H47fl (figure |21 
(e)). 

In another variant of SA-model, formula (|36|) reads as follows 



dt ] a(«y") aXgt^<g<gt^ 




(B.3) 



at < .g < g^iax ' 



In l|B.3|l . f^max'' is the maximum space gap in synchronized flow; a'^-'''™^, a!^'"°°^ ^ a^^^'^^ 
are given by (|37|l in which taken from 140(1 . 

ii:(f™o)(„^ (B.4) 

(B.5) 

where the sensitivity 



t-Csyn)/ N I -fv^ ' au (7 ^ (;£, , , 



at u < 
is:(d'='=) at w > w^, 

;;s;(frec) -g sensitivity, Vai'^"^(5) is given by A function a^^^^{g, v, vg) in (|R3)| 

is taken as follows 

a(^^)(g, t., t.,) = min(a,.ax, (F^^^) _ ^) ^ 

(B.7) 

where the function V'^^^^g) = V{g) at gma"'' < 5 < 5min°^ is a sensitivity that 

in a general case can be different from the sensitivity ^(f'''"^) in free flow. 

In the SA-model (|B.3p - l|B.7p . steady states of free flow (the curve F in figure |21 
(f)) are associated with the condition (|43|l . averaged steady states of synchronized 
flow are related to a line 5* given by the condition 

where p'^^J^^ = l/{gniax + d), and steady states for a wide moving jam are found from 
the condition ^ (figure [21 (f)). 

In contrast with the other SA-models, the SA-model (|B.3p - l|B.7p has a limited 
density range of steady states between steady states for free flow and synchronized 
flow, which are found from the condition 

a(FS)^0 ^tgt^Ug<g[^C:K (B.9) 



Eq. l(B.9p yields the following condition for these steady states in the flow-density 
plane (curve FS in flgure[21(f)) 

q^pVAp) at < p < p|:f„"\ (B.IO) 

where the density pmax'' at the maximum point for free flow is not greater than the 
density po associated with the maximum point on the curve FS (figure |21 (f)). To 
simulate a first-order F— >S transition, the steady state model solutions l|B.10|l should 
be unstable against infinitesimal non-homogeneous perturbations. This requirement 
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to the SA-model HB.3|I - HB.7() is easy satisfied through an appropriated choice of the 
function V''^^\g) and the sensitivities ^(^s)^ K'-^y'^\ In this case, numerical 

simulations of the SA-model (jB.3|) - (|B.7p made show that this model exhibits F— >S^J 
transitions in accordance with empirical results (we used the following parameters for 
the SA-model l|R3)l - ljBj)l : V{g) = Vb( tanh((5 - .go)/5i) + tanh(.go/.gi)) with Vq = 14 
m/s, go = 21 m, gi = 7 m; gtl^'^ = 24 m; A(f'^°<=) = ^(^S) ^ q.I s^i; JCCf-^c^) = 0.6 s-^; 
Xi^cc) ^ Q 4 ^(dec) jg ^^j^gj^ Tablefflwith k['^'''''^ = 1 s-^, Vc ^ 9 m/s, e ^ 
0.05; other parameters are the same as those in the SA-model (ISSt - lflT)) '). 
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